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Abstract 

g ■ Critical points that can be suppressed to zero temperature are interesting because quantum fluc- 

. tuations have been shown to dramaticahy alter electron gas properties. Here, the metal formed by 

O , Co doping the paramagnetic insulator FeS2, Fei_a;Coa;S2, is demonstrated to order ferromagneti- 



cally at X > Xc = 0.01 ± 0.005 where we observe unusual transport, magnetic, and thermodynamic 



^ ■ properties. We show that this magnetic semiconductor undergoes a percolative magnetic transi- 

o : 

\^ , tion with distinct similarities to the Griffiths phase, including singular behavior at Xc and zero 

temperature. 
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The paramagnetic (PM) to ferromagnetic (FM) transition in magnetic semiconductors is 
a prominent topic in condensed matter physics because of efforts to discover materials useful 
for spintronicsjl!]. Magnetic semiconducting materials are considered essential for use as spin 
injectors in this nascent technology, yet single phase materials with Curie temperatures, Tc, 
well above 300 K that are compatible with today's technologies have not been identified. 
Problems encountered by efforts to increase Tc of Mn doped III-V semiconductors suggest 
that a deeper understanding is needed[2| and recent theoretical investigations have provided 
some progress j^, 4]. However, disorder and strong Coulomb interactions are both central 
issues in semiconductors with magnetic impurities and this makes modeling difficult. In Mn 
doped III-V semiconductors, Mn provides localized magnetic moments as well as a smaller 
number of hole carriers that couple the local moments via the RKKY interaction. Models 
based on a magnetic-polaron Hamiltonian with random arrangements of RKKY-coupled 
moments have predicted a zero-temperature, T, percolative transition at critical magnetic 
moment and charge carrier densities^]. This critical point is governed by the competition 
between a nonmagnetic ground state and the magnetically ordered one and because of the 
doping induced disorder these materials are expected to display Griffiths phase singularities. 
Other than reports in doped LaMnOs and LaCo04 suggestive of Griffiths phase physics j^, 
there are no convincing demonstrations of Griffiths phases in magnetic semiconductors. A 
second class of materials predicted to display the Griffiths phase is the heavy fermion metals 



tuned by chemical substitution to be near a T = quantum critical point (QCP) 
In this case, rare clusters of strongly coupled magnetic moments are predicted to tunnel 
between magnetization states over classically forbidden regions resulting in non-analytic 
thermodynamic quantities at T = sl . 

In this letter we report on a magnetic and semiconducting system which resembles 
(GaMn)As with the advantage that single crystals can easily be grown and characterized; 
Co doped ton pyHte o. -fools gold'. FeS. . an .nsnlato. „.th a ^ 1 eV band gap. It i. 
isostructural to C0S2, an itinerant ferromagnet with Tc = 120 Kp^. FeS2 and C0S2 form 
a continuous solid solution over the entire concentration range, x, in Fei_2^Coa;S2, where it 



is fully spin-polarized for 0.25 < a; < 0.9 



Previous magnetic susceptibility, x, measure- 



ments found FM order at T > 2 K for x > 0.05 [10'|. Here, we discover metallic behavior 
for X > 0.001 with magnetic order a.t x > Xc = 0.01 ± 0.005 consistent with earlier work. 
In addition, we observe an evolution from a partially Kondo screened metal ai x < Xc to 
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a ferromagnet characterized by a percolative transition for x > Xc- For x ~ Xc we find a 
divergent Sommerfeld coefficient at low-T indicating non-Fermi liquid behavior along with 
a magnetic field H and T dependent magnetization M that suggests a critical transition to 
a Griffiths phase-like state. In metallic Fei_^.Co^S2 we discover evidence for the unusual co- 
existence of magnetically ordered phases and partial Kondo screening of magnetic moments. 



In contrast the model ^ for (GaMn)As predicts a Griffiths phase under the conditions of 
localized electronic carriers where magnetic-polarons are important. In light of the inferred 
competition between the RKKY coupling and Kondo screening, as well as the divergent 
thermodynamic properties at Xc and T — > 0, this system may be more closely described by 
models of /-electron materials displaying quantum criticalityjs]. 

Single crystals were synthesized by standard I2 vapor techniques from high purity starting 
materials. Crystals were etched in HCl to remove any remaining flux and characterized by 
single crystal X-ray diffraction and energy dispersive X-ray microanalysis. The Co concen- 
tration of our crystals is consistent with the saturated M at 5 T and 1.8 K and is about 70% 
of the nominal concentration of our starting materials. M and x were measured in a SQUID 
magnetometer for T > 1.8 K and a dilution refrigerator above 50 mK. The resistivity, p, 
and Hall effect were measured using four-probe lock-in techniques at 17 or 19IIz, with thin 
Pt wires attached using silver paste or silver epoxy. The specific heat was measured using a 
standard thermal relaxation method. 

Our M measurements identify magnetically ordered states with T^s shown in Fig. la 
by way of a peak in the real part of the AC susceptibility, x' (Fig. 2a), that is apparent 
in all of our samples with x > xc = 0.01 ± 0.005. We have checked that Tc signals a FM 
transition by comparing to a standard Arrott analysis of a dense set of M{H, T) data for a 
few samples with x > 0.05 fl2^. A comparison is made in Fig. la between Tc and the Weiss 
temperature, 9]y, determined from the T-dependence of x' at T's well above Tq, Fig. 2a, b. 
It is apparent that for x > Xc a. FM phase emerges and that Tc increases systematically with 
X. The scatter in these data is likely due to the variations in simultaneously grown crystals 
as is evident by the range of 6w values for crystals of the same or similar x. In Fig. lb the 
density of magnetic moments, n, determined from fits of the Curie- Weiss (CW) form (Fig. 
2b) to x{T) for T >> Tc are displayed. We have assumed an effective moment of J = 1/2 
and found n much larger than the high-if low-T saturated moment for x > x^.. This could 
indicate either a large Rhodes-Wohlfarth ratio as in itinerant ferromagnets (between 1.8 and 
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7 as compared to 3.5 in MnSi), or J > 1/2 local moments 13|. To investigate this we have 
measured the Hall effect to determine the carrier concentration, nnaii, of our samples, nnaii 
was determined at high- if to eliminate anomalous contributions. Fig. lb shows that riHaii 
ranges from 10 to 30% of x indicating that only a fraction of dopants donate electrons to 
a conducting band. Thus, it is likely that the difference between the Curie and saturated 
moments results from localized electrons with J > 1/2. 
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FIG. 1: (Color) Doping dependence (a) Tq from peak in AC susceptibility, x' from Arrott 
analysis Weiss temperature, &W: Kondo temperature, T*. Lines are linear fits, (b) Curie 
moment per formula unit from Curie-Weiss analysis and Hall carrier density per formula unit, 
UHaii, multiplied by ten. (c) Density, N, and J of spin clusters per formula unit. 

To explore more fully the properties of Fei_2,.COi.S2, we have measured the specific heat, 
C, finding nearly identical C for all our crystals above 20 K since this region is dominated by 
phonons. However, at lower T's we find a contribution that grows with x shown in Fig. 2c. 
Samples that display a finite-T peak in x also display a broad maximum in C{T)/T, albeit 
at a lower T. The idea of spin clusters suggested by M{H,T) and x'{T) was probed by 
comparing C{T) and x'{T) to determine J, the fluctuating moment above Tq. At T > 2K 
we fit C{T) by the sum of 7T, a phonon term jST^ {(3 = 2.17 x 10'^ J/mol K^), and a 
Schottky term due to localized magnetic moments with nJ{J + 1) determined by fits of the 
CW form to x'{T) (Fig. 2c). The best fit n and J are shown in Fig. Ic where we observe 
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FIG. 2: (Color) Susceptibility and Specific heat, (a) Temperature T dependence of real part of AC 
susceptibility, x' ■ Inset: x'{H) at 4 K with power law fit (red line), (b) T dependence of 1/x' at 
and 1 mT. Red line is fit of Curie- Weiss form at high-T. Inset: l/x'(0) vs. reduced T, (T/Tq — 1). 
Red line is 1/x' = a{T-T^)^ with 6 = 0.45 ±0.04 and = 24 ±0.5 K. (c) Specific heat, C divided 
by T, vs. T for X = 0.002 (circles), 0.005 (bullets), 0.007 at if = (blue squares), H = IT (solid 
line), and ii = 3 T (dashed line), 0.03 (FM) (triangles), and 0.045 (FM) (diamonds). Red line is 
fit of the form aT~" with a = 0.69 ± 0.05, for x = 0.005. Dotted lines are fits of a Sommerfeld 
plus Schottky model to the data at T > 2 K. 

that J > 1/2, n < X, and that both grow with x. Thus, spin clusters are consistent with 
both C{T) above 2 K and x'{T). 

The magnetic transitions were explored in detail by measuring the low- if x'(if, T) as in 
Fig. 2a, b. What is interesting is that both the T and magnitude of the x' niaxima are 
significantly suppressed by very small H. x'{H) at 4 K is displayed in the inset to Fig. 2a 
where a power law form x'{H) = hH~^ with (3 = 0.62 ± 0.03 describes the data well. In 
addition to this extreme field sensitivity we have observed deviations from CW behavior. In 
Fig. 2b we plot 1/x' at DC fields of and 1 mT with dramatic changes evident below 30 
K. The 1 mT data follow a CW form, x' oc (T - OwV^ with 9w = IS K for T > Ow In 
contrast, the H = data cannot be described by a CW form to much higher Ts and the 
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deviation is in the direction of smaller 1/x'- The growth of x' beyond the CW form is an 
indication of short range FM correlations [See e.g. [91]] as spin clusters imply larger x'- In the 
inset of Fig. 2b we display a fit of the form 1/x' = (T — TqY to the data, Tq = 24 ±0.5 
K ^ 2Q\Y the critical T for the largest clusters and 6 = 0.45 ± 0.04, suggestive of Griffiths 
phase formation. 

Samples that remain PM down to our lowest T have a C{T)/T that decreases with T 
and that can be accurately described by a aT~" form with a < 1. The red line in Fig. 
2c represents this form with a = 0.69 ± 0.03. Application of if > suppresses the low-T 
C{T) /T of all of our samples and the PM samples display a C{T)/T that resembles our FM 
samples. If we make the assumption that the conduction electron gas acts independently 
from a set of weakly interacting local moments, then C{T) would be well fit by the sum of 
7T, j3T^ and a Schottky-like anomaly down to zero-T. At finite field, the anomaly would 
evolve into a Schottky peak. Although this description works well above 1 K, our low- 
T data does not conform to this simple picture. Interestingly, materials in proximity to 
metal-insulator transitions have a diverging C{T)/T also described by a T~" form. This is 
ascribed to the random position of local moments that are interacting antiferromagnetically 
leading to a singlet ground state However, since our samples are either FM or nearly 
FM, we do not consider this to be a likely explanation of our data. Instead, we suggest that 
the divergence of C{T)/T for a; ~ Xc is indicative of Griffiths phase physics and/or to the 
proximity to a FM QCP. Magnetic fields or ordering return C{T)/T to a Fermi liquid form 
at T's proportional to if or Tc. 

The charge carrier transport properties of our crystals are presented in Fig. 3. While 
nominally pure FeS2 displays insulating behavior, Co doping at a level of x = 0.001 is 
sufficient to create a metal (Fig. 3a). Larger x tends to increase nnaii and decrease the 
resistivity, p. One interesting feature is that the T-dependence of p for x ~ Xc closely 
follows a form with ^ = 1.6 ±0.1 over a wide T-range similar to that found in MnSi near 
the critical pressure for suppression of ferromagnetism. It is also consistent with the spin 
fiuctuation model of nearly FM metals where a T^/^ dependence is predicted for itinerant 
magnets. However, as noted above, the magnetic moments are likely to result from localized 
electrons so it is not clear that this model is appropriate. Evidence for the importance 



of magnetic fiuctuations 
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16| in determining p(T) can be seen in Fig. 3b where the 



T-derivative of p normalized by po = p(4i^), dp/dT/po, is displayed. The normalization 
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removes error in crystal geometry determination as well as changes due to variations in 
nHaii- This quantity is strongly peaked near Xc- Thus, we observe a reduced power-law 
behavior of p{T) along with an increased scattering rate over a wide T-range in proximity 
to the zero-T critical point for magnetism. 
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FIG. 3: (Color) Carrier Transport (a) Resistivity, p, vs. temperature to the 1.6 power, T^'^, for 
a subset of our crystals, (b) T-derivative of p, dp/dT, normalized by p(4 K), pQ, vs T. (c) T 
dependence of p for x = 0.005 at magnetic fields, H, identified in the figure. Blue line is fit of a 
In T behavior for H = and red line is fit of Kondo theory [f?!. (d) Scaling plot of magnetoresistance 
at Ts and Hs identified in (c) and (d). 



In addition, a second contribution to p is apparent at T < 20 K as demonstrated for a 
PM sample in Fig. 3c. Here p increases with decreasing T in a manner that is well described 
by a logarithm over more than a decade in T with Kondo temperatures T* = 0.8, 1.5, 1.4, 
and 2.5 K for x =0.002, 0.004, 0.005, and 0.007, Fig. lapj. We have also measured a 
large negative magnetoresistance (MR) (Fig. 3c, d) that is identical in the transverse and 
longitudinal current directions. This indicates a spin, rather than orbital, mechanism for 
the MR consistent with a Kondo effect dominating p(T, H). Furthermore, all of our T- and 
if-dependent data can be scaled by a single ion Kondo form; p(T, H) — piT, 0)/p(T, 0) scales 
as H/ (T + T*) for T > T* as shown in Fig. 3djl8|. We conclude that a single energy scale, 
likely a Kondo coupling of conduction electrons with the local moments associated with the 
Co ions, determines the low-T p of Fei_xCo^S2. 

The observation of power-law divergent C and x' for a; ~ Xc with similar exponents 
suggests that a single physical mechanism describes both. The sensitivity to magnetic fields 
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indicates that there is weak couphng between clusters creating a fragile magnetic order for 
small X. In Fig. 4 we summarize the evolution of Fei_2:Coa;S2 from a strongly PM low- 
carrier-density metal to a FM with increasing x. In Fig. 4a, for x < Xc and T = 0, the 
system consists of a small density of magnetic moments localized on the Co impurity sites. 
The electron carriers screen only a portion of these moments via the Kondo coupling evident 
in the carrier transport. As x increases to Xc in Fig. 4b there is a percolative transition at 
T = and we observe a divergent C{T)/T at low-T. Clearly for a magnetic transition to 
occur at X ~ 0.01 long range interactions between moments is necessary. Thus, a percolative 
transition occurs at Xc << 0.2 required for a face-centered-cubic lattice with short range 



interactions [19|]. At still larger x. Figs. 4c and d, the FM state is complete at low-T. For 
Ts above the critical point the system is dominated by super-paramagnetic regions where 
magnetic moments are correlated, but long range ordering is not apparent. The system 
consists of weakly coupled clusters that are easily modified by small H. In addition, the 
clusters have a finite probability to tunnel to nearly degenerate M states. 
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FIG. 4: (Color) Evolution, (a) Low Co concentration (x < Xc) with disperse local moments (red 
arrows) and low density of itinerant electrons (blue arrows). Electrons partially screen moments 
via Kondo coupling, (b) At percolative magnetic transition occurs at T = 0. Ordered 

regions are red with magnetization, M, direction green, (c) Larger x, x > Xc, fully ordered at 
T = with large M domains, (d) T > Tc, clusters of strongly coupled spins form for T < T^. 
Tunneling of clusters indicated by the double-ended green arrows. 
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where disorder is sufficient 



These features are described by Griffiths phase models 
to cause clusters of localized short lived magnetic order at T's above the global ordering 
temperature. As the system is cooled toward Tc, these clusters grow and display switching 
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of M via tunneling. The consequence of this model is a power-law form, 1/(T — T^Y~^, 
with A < 1, of the thermodynamic quantities above Tc that can be suppressed with small to 
moderate magnetic fields. For Fei_j.Co^S2 0.30 < A < 0.55. Griffiths phase physics has been 
suggested to explain the non-Fermi liquid behavior of heavy fermion antiferromagnets with 
Neel temperatures driven to zero by chemical substitution. However, Millis et al. 20|] point 
out that the coupling of magnetic clusters to conduction electrons will suppress tunneling, 
removing the non-Fermi liquid response in the model. In contrast to well developed con- 
ductors, our materials are nascent metals with poor electrical screening formed by doping 
an insulator. Although Fei_sCo2,.S2 does not appear to be described by a magnetic-polaron 
Hamiltonian as in Ref. [3], our data suggest that Griffiths phase anomalies, likely influenced 
by Kondo screening, are in fact observable in magnetic semiconductors. The result is singular 
behavior at T = similar to models of Griffiths phases in /-electron antiferromagnets. 

We thank I. Vekhter and C. Capan for discussions. JFD, DPY, and JYC acknowledge 
support of the NSF under DMR0406140, DMR0449022, and DMR0237664. 
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